Search for a pentaquark decaying to pK 
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Abstract 

We present a search for a pentaquark decaying strongly to pK$ in ^yN collisions at 
a center-of-mass energy up to 25 GeV/c 2 . Finding no evidence for such a state in 
the mass range of 1470 MeV/c 2 to 2200 MeV/c 2 , we set limits on the yield and on 
the cross section times branching ratio relative to S*(1385) ± and K*(892) + . 
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1 Introduction 

Nearly 30 years ago Jaffe proposed the existence of bound (mass below thresh- 
old for strong decay) multiquark states including QQqq states and the H di- 
hyperon [1] based on calculations using the bag model [2] . As the years passed 
and no convincing evidence for non mesonic and non baryonic states was found 
the field languished. 
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Between January 2003 and March 2004, however, the pentaquark field was 
reenergized when no less than ten independent pentaquark observations at a 
mass around 1540 MeV/c 2 were reported [3,4,5,6,7,8,9,10,11,12]. The presumed 
quark content of the reported states was (suudd). The reported widths were 
consistent with detector resolution even though a strong decay is allowed. 
Even more amazing was the prediction six years before by Diakonov et al. 
of just such a state at a mass around 1530 MeV/c 2 and a width less than 
15 MeV/c 2 [13]. Observations of a doubly-strange pentaquark [14] and a charm 
pentaquark [15] have also been reported by single experiments. 

The original observations were made in the nK + mode. For a state to decay 
strongly to nK + , it must be composed of at least 5 quarks. Other observations 
have been made in the pK s mode which is not manifestly exotic since the K s 
can originate from a K° or K°; a pK° decay is not exotic while a pK° decay is 
exotic. Since that time, many other experiments have failed to find evidence of 
pentaquarks [16,17,18,19,20,21,22,23,24,25,26]. Most of these experiments are 
higher statistics and higher energy than the observing experiments and gener- 
ally search the pK s decay mode. Recently CLAS, which previously reported 
two observations failed to find pentaquarks in a third attempt [27]. This letter 
presents a search of the FOCUS data for the + (suudd) pentaquark candi- 
date in the decay mode + — > pK s . Cross sections will be measured relative 
to three well known states with a similar decay topology: S*(1385) ± — > 
(two states) and K*(892) + -> K° s n + 2 . 



2 Event reconstruction and selection 

The FOCUS experiment recorded data during the 1996-7 fixed-target run at 
Fermilab. A photon beam obtained from bremsstrahlung of 300 GeV elec- 
trons and positrons impinged on a set of BeO targets. Four sets of silicon 
strip detectors, each with three views, were located downstream of the tar- 
gets for vertexing and track finding. For most of the run, two pairs of silicon 
strips were also interleaved with the target segments for more precise ver- 
texing [28]. Charged particles were tracked and momentum analyzed as they 
passed through one or two dipole magnets and three to five sets of multi- 
wire proportional chambers with four views each. Three multicell threshold 
Cerenkov counters, two electromagnetic calorimeters, and two muon detectors 
provided particle identification. A trigger which required, among other things, 
>25 GeV of hadronic energy passed 6 billion events for reconstruction. 

The data used for this analysis come from a subset of FOCUS data which con- 
tain vee candidates (Kg — > it + it~ and A — > pn~). There are four vee candidate 

2 Charged conjugate states are implied unless explicitly stated otherwise 
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types which are used in this analysis. SSD vees decay in the vertex region and 
have decay tracks which are found in the silicon system. Magnet vees decay 
further downstream and tracks are reconstructed in the wire chambers. The 
magnet vees are divided into three types SS, TS, TT for stub-stub, track-stub, 
and track-track depending on whether the decay particles are tracked in only 
the upstream 3 wire chambers (stubs) or in all 5 wire chambers (tracks). A 
full description of vee reconstruction in FOCUS can be found in Ref. [29]. 

The vee selection requires a reconstructed vee vertex with some quality, parti- 
cle identification, and mass cuts. The quality requirements include good track 
quality, a good vertex for the two tracks, and a minimum vee momentum 
of 5 GeV/c. The mass requirement for a Kg candidate is to have a normal- 
ized mass within 4 (5) er of the nominal mass for SSD (magnet) vees. The 
A mass requirement is for the invariant mass to be between 1.09 GeV/c 2 
and 1.14 GeV/c 2 . The particle identification cuts on the vee daughters use 
the FOCUS Cerenkov identification algorithm [30]. This algorithm returns 
negative log-likelihood (times two) values Wi(j) for track j and hypothe- 
sis i G {e, 7i, K, p} based on the light yields in the phototubes covering the 
Cerenkov cone of the track. The information from all three Cerenkov detectors 
is combined. The vee daughter pion candidates must not be strongly incon- 
sistent with the pion hypothesis: W m i n (7r) — W vr (7r) > —5 where W m i n is the 
minimum of the four hypotheses. The small phase space of the A — > pn~ 
decay and the forward nature of the FOCUS spectrometer require the proton 
to have a higher lab momentum than the pion. Thus, the higher momentum 
track is chosen to be the proton and it must have W min (p) — W p (p) > —5 
and W 7r (p) — yV p (p) > 1. To reduce combinatorics, the event is only kept if 
the number of vees passing the above cuts was no more than two. At least 
one and no more than seven good quality charged tracks with momentum 
greater than 5 GeV/c must be found in addition to the vee(s). These tracks 
and any SSD vees must be consistent with originating from a single vertex 
with confidence level greater than 1%. The total number of good quality vees 
plus good quality charged tracks must exceed two. After applying all of the 
above cuts, the vee sample within 3a of the nominal mass is selected which 
contains 72 million K s — > 7r + 7r~ (9 million A — > ptt~) candidates of which 
90% (95%) are signal as shown in Fig. 1. Each vee in an event is combined with 
the good quality charged tracks in the event to search for K*(892) + — > K s 7r + , 
S*(1385) ± — > A ^, and + — > pK s . The Cerenkov requirements on the 
charged tracks are W m i n (7r) — W^iix) > —2 for pions plus Wk{^) — W^tt) > 1 
for the pion from the S*(1385) ± decay. The proton from the + decay must 
have W n (p) - W p (p) > 8 and W K (p) - W p (p) > 3. The Cerenkov cuts were 
optimized using signal Monte Carlo with data background. The proton re- 
quirements, in particular, are very stringent and reduce the misidentification 
rate to nearly zero. 

The K s ir + and A ^ are shown in Figs. 2 and 3, respectively. In both cases, the 
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Yield: 67877952 ± 11430 
f S/B (±3a) = 8.9 




[ M(7t + 7t~) - M(K S ) ] / CT, 



-3-2-1 i ; 
[ M(pTt-) - M(A) ] / a K 



Fig. 1. The normalized mass plots of K$ — ► tt^tt and A 
Events inside the vertical lines are selected for analysis. 
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Fig. 2. K*(892) + fit with an S-wave Breit-Wigner and combinatorial background. 

signal was best fit with an S-wave Breit-Wigner with an energy independent 
width even though a P-wave energy dependent width would be more appro- 
priate. The Breit-Wigner was convoluted with a Gaussian for the detector res- 
olution obtained from a Monte Carlo simulation. The i^*(892) + , S*(1385) + , 
and E*(1385)- resolutions are 5.1 MeV/c 2 , 3.2 MeV/c 2 , and 3.2 MeV/c 2 . 
The background was fit to the form aq b exp (cq + dq 2 + eq 3 + fq A ) where a — / 
are free parameters and q is the Q-value (invariant mass minus component 
masses) . 



3 Pentaquark search results 



The pKg and pKg invariant masses are plotted using the standard selection 
criteria in Fig. 4. There are no significant differences between the two charge 
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Fig. 3. X*(1385) fits with an S-wave Breit-Wigner and combinatorial background. 

states so for the remainder of the analysis we combine the charge conjugate 
states. The combined sample with standard cuts and with an additional mo- 
mentum asymmetry cut is plotted in Fig. 5. The momentum asymmetry cut, 
requiring the proton to have a higher momentum than the K s in the pen- 
taquark decay, has been suggested as a method of reducing background. The 
true effect of the cut is to sculpt the mass distribution into a more peaked 
structure near the location of the previously observed pentaquarks and is not 
used in the analysis. In Fig. 6 the total sample is fit to a background curve of 
the form aq b exp (cq + dq 2 + eg 3 + fq A ) where a — / are free parameters and q 
is the Q- value: q = M(pK$) — m p — m^o. No evidence for a pentaquark near 
1540 MeV/c 2 or at any mass less than 2400 MeV/c 2 is observed. To set a limit 
on the yield we need to make some assumptions about the width of the state. 
We consider two cases: one with a natural width of and one with a natural 
width of 15 MeV/c 2 . In the first case, the signal is fit with a Gaussian with a 
width set by the experimental resolution. In the second case, the signal is fit 
with an S-wave Breit-Wigner with an energy independent width convoluted 
with the experimental resolution. The experimental resolution in MeV/c 2 is 
approximately <r(MeV/c 2 ) = —22.5 + 19.48m — 1.99m 2 where m is the mass 
in GeV/c 2 . 

A series of 731 fits to the observed pK$ mass plot were performed using 
the background and signal shapes described above. The signal mass is varied 
in 1 Me V/c 2 steps from 1470 to 2200 MeV/c 2 and a binned log-likelihood fit 
using Minuit [31] is performed. The ±1<t errors are defined as the point where 
Alog£ = 0.50 relative to the maximum log£, while continually adjusting the 
background parameters to maximize log C The 95% CL lower limit is defined 
similarly with Alog£ = 1.92. Both are obtained using Minos [31]. The 95% 
CL upper limit is constructed as follows: The likelihood function C versus yield 
is determined by maximizing log£ for many different (fixed) yields, allowing 
background parameters to float. The likelihood function is integrated from a 
yield of to oo to obtain the total likelihood. The 95% CL upper limit on the 
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Fig. 4. Invariant mass distribution of pKg separated by charge. Standard cuts are 
applied. 
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Fig. 5. Invariant mass distribution of pK$ for both charge states. Solid histogram 
shows the result for standard cuts and the dashed histogram is with an additional 
cut requiring the proton momentum be greater than the K$ momentum. 

yield is defined as the point where 95% of the total likelihood is between a yield 
of and the upper limit. This definition of an upper limit is used rather than a 
counting based Feldman-Cousins type limit due to the large background which 
results in Gaussian errors. The fitted yield, 1-a errors, and 95% CL limits are 
shown in Fig. 7. Of the 1462 fits, none of them finds a positive excursion 
greater than 5a. Previous pentaquark observations have occurred between 
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Fig. 6. Invariant mass distribution of pKa for both charge states with standard cuts. 



1520 MeV/c 2 and 1555 MeV/c 2 with a very small natural width. In this region 
the largest excess we find is a 2.5a excess at 1545 MeV/c 2 for the T = MeV/c 2 
fits. A 3.1cr excess is seen at the same location for the T = 15 MeV/c 2 but this 
width would be inconsistent with the previous observations and is certainly not 
a convincing observation. Additionally, this excess occurs in a region where 
the background distribution is peaking which can give the appearance of a 
signal. Given the large number of fits, the appearance of 2-3(7 excesses is not 
unlikely. 



To compare with other experiments, the limits on yield must be converted 
to limits on production times (unknown) branching ratio. We choose to nor- 
malize the 6 + production cross section to S*(1385) ± and i\"*(892) + because 
the reconstructed decay modes of these particles E*(1385) ± — > A ^ and 
K*(892) + — * Kgir + are very similar, in terms of topology and energy release, 
to the signal. Thus, we attempt to determine 



a (0 + ) ■ BR (9" 



>pK§) 



a (if* (892)- 



and 



a(0+) ■ BR (Q+->pK%) 
a(E*(1385)±) 



Rewriting Eq. 1 in terms of measured yields (Y) and efficiencies (e), we find: 
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Fig. 7. Pentaquark yields and upper limits. Top (bottom) plots show results for a 
natural width of (15 MeV/c 2 ). The shaded region includes the 1 a errors with the 
central value in the middle. The outer curves show the upper and lower limits. 

a(6 + )-BR(0+^^) Y(G+)-BR(G+^pK s )-e K * {892)+ 



a (if* (892)+) 

a (9 + ) • BR (6 + ^pKg) 
a(S*(1385) ± ) 



e e+ ^ K o -Y(K*(892)+) 

Y(e + )-BR(9 + ->^g)-6 S »( 1385 ) ± 
e e+ ^o -Y(S*(1385)±) 



(2) 



All of the efficiencies include the reconstruction and selection efficiencies plus 
corrections for unseen decays of parent particles. The G + — > pK$ efficiency 
only includes the correction for the unseen Kg decays, not corrections for 
©+ — ► pKl or other G + decay modes such as + — > nK + . The K*(892) + and 
S*(1385) efficiencies include all branching ratio corrections from the PDG [32]. 
For K*(892) + this corrections comes from BR(K S -> tt+tt-) = 0.6861, BR(K° - 
K° s ) = 0.5, and BR(lf(892) + -> K°tt + ) = 0.667 while for S*(1385) ± the rel- 
evant branching ratios are BR(A° — > p7r~) = 0.64 and BR(S*(1385) ± — > 
A°7T ± ) = 0.88. Determining reconstruction and selection efficiency (including 
acceptance) is described below. 

The FOCUS detector is a forward spectrometer and therefore acceptance de- 
pends on the produced particle momentum. The production characteristics 
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of the pentaquark are the largest sources of systematic uncertainty in this 
analysis. We choose a particular production model to obtain limits and pro- 
vide sufficient information about the experiment for other interested parties 
to obtain limits based on other production models. The production simula- 
tion begins with a library of e~ and e + tracks obtained from a TURTLE 
simulation [33] of the Wideband beam line. From this library, an individual 
track is drawn and bremsstrahlung photons created by passage through a 
20% Xq lead radiator. Photons with energy above 15 GeV are passed to the 
Pythia [34] Monte Carlo simulation. The Pythia simulation is run using 
minimum bias events (MSEL=2) with varying energies (MSTP(171)=l). Options 
controlling parton distributions and gluon fragmentation were set to avoid 
heavy quark production (MSTP(58)=3 and MDME(156~ 160, 1)=0). Also the 
center-of-mass minimum energy cut off was reduced from 10 GeV to 3 GeV 
(PARP(2)=3). However, the minimum photon energy requirement gives an ef- 
fective minimum center-of-mass energy of 5.3 GeV. Since Pythia does not 
produce pentaquarks, another particle must be chosen to represent the pen- 
taquark. Other than mass, the most important effect on the production is 
the number of quarks a particle has in common with the initially interacting 
hadrons, due to the nature of the Pythia string fragmentation model. The 
H*(1530)° and S*(1385) + particles are chosen to represent the extremes in 
the production of a pentaquark. The E*°(ssu) (£ + (swu)) can obtain at most 
33% (67%) of the remaining quarks from the target nucleon valence quarks, 
while the Q + (suudd) can take 60%. In all cases, the charge conjugate particles 
must obtain all quarks from the vacuum. The mass of the particle chosen to 
represent the pentaquark, S*(1530)° or S*(1385) + , is set to the appropriate 
value in PYTHIA, by setting PMAS( 190,1) or PMAS (187,1), respectively. 

To calculate the relative cross sections in Eq. 2 we need efficiencies for S*(1385) ± — > 
A°7T ± , K*(892) + -> K° s 7r + , and 0+ -> pK° s . These efficiencies are obtained 
from the FOCUS Monte Carlo simulation. The dominant uncertainty in the 
efficiency determination is the modeling of the production characteristics of 
the parent particle. For the observed particles, S*(1385) ± and K*(892) + , we 
can compare the data and Monte Carlo directly and adjust the Monte Carlo 
simulation to produce the correct data distribution. Even this is not suffi- 
cient, however, because areas where the efficiency is zero cannot be accounted 
for. For S*(1385) ± and fT*(892) + , we run a weighted Monte Carlo simula- 
tion which matches the Monte Carlo momentum distribution with the ob- 
served data momentum distribution in the region for which the acceptance is 
not zero. The dominant source of uncertainty for the S*(1385) =t — > Air^ 1 and 
K*(892) + — > Kgn + efficiencies is the our lack of knowledge of the fraction of 
events completely outside of our acceptance (momentum less than 15 GeV/c). 
The weighted Pythia Monte Carlo predicts 67% (79%) of the A"* (892)+ 
(S*(1385) ± ) particles are produced with momentum less than 15 GeV/c. To 
obtain an estimate of the efficiency uncertainty, we assume that the number 
of particles with momentum less than 15 GeV/c can be off by up to a factor 
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Fig. 8. Acceptance versus mass for pentaquark candidates. Upper (lower) curve is for 
a pentaquark produced as a H*(1530)° (£*(1385) + ) and the middle is the average. 

of 2 (high or low). This leads to a relative uncertainty on the K*(892) + and 
S*(1385) ± efficiency of 43% and 49%, respectively. The G + — > pK$ efficiency 
is taken as the average of the efficiencies obtained from using £*(1385) + and 
H*(1530)° as the substitute particle while the uncertainty is half the difference 
between the two efficiencies. The @ + — > pK$ with Kg — ► 7r + 7r~ efficiency 
versus mass (with no branching ratio corrections) is shown in Fig. 8. The av- 
erage uncertainty in e e +^ pK o is approximately 26%. It may seem incongruous 
that the relative uncertainty of the efficiency of an unknown particle (~26%) 
is less than that for the high statistics normalizing modes (>40%). The effi- 
ciency uncertainty of the high statistics modes reflects the lack of knowledge 
of production outside of our acceptance. However, it is reasonable to assume 
that discrepancies in the Monte Carlo simulation will be similar for the signal 
mode and the normalizing mode and therefore adding the uncertainty to the 
signal mode is double-counting. Note that the signal and normalizing efficien- 
cies only appear as a ratio. 

We also report the relative cross sections in the region where our acceptance is 
good, that is for parent particle momenta greater than 25 GeV/c. This dramat- 
ically reduces the systematic uncertainties associated with the measurement. 
The uncertainty due to the production of S*(1385) ± and K*(892)° is minimal. 
The uncertainty in the + efficiency is also dramatically reduced from approx- 
imately 26% to about 6% as shown in Fig. 9. The number of reconstructed 
fT*(892) , £*(1385)+, and E*(1385)~ at momenta greater than 25 GeV/c is 
7.88 million, 127,000, and 212,000, respectively compared to the total sample 
of 8.98 million, 151,000, and 256,000, respectively with no momentum cut. 

The upper limit on the yield was obtained by mathematically integrating the 
likelihood function from to oo and then integrating from to 95% of the 
total likelihood integral gave the 95% CL upper limit. To obtain the limit on 
cross section requires a different approach due to the significant systematic 
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Fig. 9. Acceptance versus mass for pentaquark candidates. Lower (upper) curve 
is for a pentaquark produced as a H*(1530)° (S*(1385) + ) and the middle is the 
average. The pentaquark is produced and reconstructed with momentum greater 
than 25 GeV/c. 



uncertainties. We use a method based on a note by Convery [35] which is 
inspired by the Cousins and Highland [36] philosophy for including systematic 
uncertainties. The Cousin and Highland prescription is appropriate for low 
background experiments with Poisson errors while the Convery proposal is 
applicable to the Gaussian errors which result from the large background in 
our case. Modifications to the Convery approach are made to give an exact 
solution [37]. 

Before systematics are considered, an analysis using a maximum likelihood fit 
returns a central value for the branching ratio (B) and a statistical error (<Jb)- 
The likelihood function is 



p(B) oc exp 



[B-B) 
2al 



(3) 



Following the notation in Convery, we associate S with the nominal efficiency 
and ag as the error on the efficiency Adding the uncertainty on the efficiency 
changes the likelihood to: 



p(B) oc / exp 
Jo 



-{SB/S-Bf 
2^1 





'-{S-Sf 


exp 


1*1 \ 



dS 



(4) 



Using Mathematica® , removing unimportant multiplicative constants, and 
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changing variables from as to a e = a s /S, the integral in Eq. 4 becomes: 



p(B) oc 



B2 , 1 

-3 — r —jI 



cxp 



— (B — B) 2 
2{B*a* + o%) 



erf 



erf 



BBa 2 + a 2 



V2a e a B ^fB 2 a 2 + a 2 B _ 



l)a 2 B - Ba 2 (B - BS) 



V2Sa e a B JB 2 a 2 + a 2 B 



(5) 



We integrate Eq. 5 from to 00 to obtain the total probability and then 
integrate from to the point at which 95% of the total probability is included 
and define this as our 95% CL upper limit. The branching ratio B of Eq. 5 
is simply the relative cross sections times the unknown pentaquark branching 
ratio as in Eq. 2. The relative uncertainties on the efficiency for the signal and 
normalizing mode are added in quadrature to become a e in Eq. 5. Furthermore, 
S is the relative efficiency between the signal and normalizing modes and <jb is 
the statistical uncertainty on the branching ratio due simply to the uncertainty 
in the signal yield. 



Figure 10 shows the results for J(k^(^2)+) K S ^ w ^ n an assumed natural 
width of (15) MeV/c 2 for the top (bottom) plot. This is the result corrected 
for all undetected particles. The shaded band shows the ±la limits with statis- 
tical uncertainties only; the line in the middle of the band is the central value. 
The top curve shows the 95% CL upper limit using the method described 
above including statistical and systematic uncertainties. The curve between 
the full upper limit and the la band is the 95% CL upper limit using the 
method described above with no systematic uncertainties included. The large 
systematic uncertainties are due to the attempt to correct for the vast major- 
ity of particles outside of our acceptance. While this systematic uncertainty 
significantly degrades the limit, the production times branching ratio of the 
pentaquark relative to K*(892) + production is still less than 0.0013 (0.0032) 
(95% CL) over the entire mass range for a natural width of (15) MeV/c 2 . 
Figure 11 gives the same results for ~^ K s) _ j n cage ^ 95% 

CL upper limit on the pentaquark production times branching ratio relative 
to S*(1385) ± is 0.025 (0.062) over the entire mass range for a natural width 
of (15) MeV/c 2 . 



Figures 12 and 13 show the same results for the restricted range of momentum 
greater than 25 GeV/c. That is, they show limits on relative cross sections for 
particles (9+ 7^(892)+, S*(1385) ± ) produced with p > 25 GeV/c. 
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Fig. 10. 2 (^^(892 ) + versus mass. Top (bottom) plots show results for a + 
natural width of (15 MeV/c 2 ). The shaded region encompasses the 1 a statistical 
uncertainty with the central value in the middle. The top curve shows the 95% CL 
upper limit including systematic uncertainties while the middle curve is the 95% 
CL upper limit with statistical uncertainties only. 

4 Conclusions 



We find no evidence for pentaquarks decaying to pK$ in the mass range of 
1470 MeV/c 2 to 2200 MeV/c 2 . In contrast, we observe 9 million K*(892) + 
K^7T + particles and 0.4 million S*(1385) ± — > A ^ particles which have a 
very similar topology and energy release. We set 95% CL upper limits on the 
yield over the entire mass range with a maximum of 1300 (3000) events for an 
assumed natural width of (15) MeV/c 2 . We also obtain 95% CL upper limits 
on the cross section for pentaquark production times the branching ratio to 
pK° s relative to 7C(892)+ K° s n+ and £*(1385) ± -> A vr ± . These limits 
are determined for two cases. The first case is for parent particles produced 
at any momenta (albeit with a minimum center-of-mass energy of 5.3 GeV) 

where we find a maximum upper limit of }(K^(m2)+) K ^ < 0-0013 (0.0033) 

and "^l^lss^^ < 0- 023 (0-057) at 95% CL for a natural width of 
(15) MeV/c 2 . In the second case we measure the relative cross sections for 
parent particles with momenta above 25 GeV/c (a region of good acceptance) 

and calculate 95% CL limits of ^w^Sfo < 0.00012 (0.00029) and 
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Fig. 11. CT (s*^3g^+)^^(x:* (^85)-) versus m ass. Top (bottom) plots show results for 
a + natural width of (15 MeV/c 2 ). The shaded region encompasses the 1 a 
statistical uncertainty with the central value in the middle. The top curve shows 
the 95% CL upper limit including systematic uncertainties while the middle curve 
is the 95% CL upper limit with statistical uncertainties only. 

^ff^S^^ < 0.0042 (0.0099) for a natural width of (15) MeV/c 2 . 

Very few of the observing experiments report results for A^*(892) + — > KgTT + 
and E*(1385) ± — > A ^ yields. One CLAS result apparently finds a yield of 
~1000 A"* (892)+ [38] while SVD reconstructs ~125 K*(892) + -> K° s ir+ and 
-100 £*(1385)+ -> A°7T + decays [10]. The FOCUS results presented here 
represent samples that are more than 2000 times larger. Unfortunately, dif- 
ferences in production between (mostly) low energy experiments which have 
reported observations and a high energy experiment such as FOCUS prevent 
any definitive conclusions from being drawn. 
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